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Transcriptional repression may cause transcriptional noise by a competition between repressor and 
RNA polymerase binding. Although promoter activity is often governed by a single limiting step, we 
argue here that the size of the noise strongly depends on whether this step is the initial equilibrium 
binding or one of the subsequent uni-directional steps. Overall, we show that non-equilibrium 
steps of transcription initiation systematically increase the cell to cell heterogeneity in bacterial 
populations. In particular, this allows also weak promoters to give substantial transcriptional noise. 


I. INTRODUCTION 

Protein production in living cells is the result of 
the combined dynamics of transcription and translation, 
through the activity of first RNA polymerase (RNAP) 
that synthesizes messenger RNA (mRNA), and subse¬ 
quently the ribosomes that translate the information on 
mRNA to proteins. Because each mRNA typically are 
translated many times [l], the fluctuations in protein 
number are sensitive to fluctuations in the number of 
produced mRNA Si- Therefore there have been sub¬ 
stantial interest in determining the noise in this number 
EH3- This noise is primarily governed by the stochastic 
dynamics of RNAP around the promoters, which are the 
regions on the DNA that direct initiation of the tran¬ 
scription process. 

With recent availability of the technology for counting 
individual mRNAs in E. coli cells EH3 , it has become 
feasible to quantify the interplay between noise in gene 
expression and dynamics around the promoter. The de¬ 
gree of cell-to-cell variability in the number of a given 
mRNA is often quantified by the Fano factor, the ra¬ 
tio between the variance and the mean. The Fano fac¬ 
tor exceeds one when the transcription is bursty. Such 
transcription burstiness can be obtained from a model 
where a gene switches between an “on-state” with high 
promoter activity and an “off-state” with low activity 
Ellil- In this simple scenario, Fano factor gives the 
estimate of the number of transcripts produced per “on- 
state” , and such a scenario can be realized by different 
molecular mechanisms. 

Transcriptional regulators influence RNA polymerase 
(RNAP) access to promoters, and may cause alternat¬ 
ing periods of low and high promoter activity, depend¬ 
ing on the presence or absence of the regulator near the 
promoter (Fig. [TJ. When the repressor is the source of 
the burstiness, the measurements of the Fano factor for 
mRNA levels may allow for quantification of the relative 
sizes of on-rates of transcriptional repressors and on-rates 
of RNAP [sj. A recent study (To| reported Fano factors 
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FIG. 1. Illustration of the competition between a simple tran¬ 
scriptional repressor (blue) and the RNA polymerase (red) in 
terms of the time intervals they occlude the promoter. No¬ 
tice that a bound RNAP takes time to initiate transcription. 
It is only when the promoter is open, that there is a direct 
competition for the available space. The probability that the 
repressor wins this competition is kb/{kb + r on ) where r on is 
the effective on-rate of the RNAP (see Fig. 2) and kb is the 
binding rate of the repressor. The number of times the RNAP 
binds before the repressor rebinds is given by r on /kb. 


in the presence of a transcriptional repressor. The mea¬ 
sured dependence of noise on repressor concentration was 
reproduced by using a one-step model for transcription 
initiation, assuming that RNAP binding to the promoter 
sequence is the rate limiting step. 

However, transcription initiation in Escherichia coli in¬ 
volves at least three steps: closed complex, open complex 
and elongation initiation 00 (see Fig. [21 of which the 
two later steps are often limiting 0 , 0 - fl5| . Measur¬ 
ing the distributions of time intervals between two sub¬ 
sequent transcription events directly demonstrated that 
the tetA promoter has at least two limiting steps 0]. 
In cases where promoter activity is limited by later steps 
of the initiation process, the RNAP is bound to the pro¬ 
moter for a longer period. This inhibits the access for 
subsequent RNAPs as well as for transcription factors in 
the occluded region 01 as indicated by the red squares 
in Fig. |T] In fact ref. 0 studied a synthetic model sys¬ 
tem where the time RNAP spends on a promoter allowed 
a four fold repression of a partly overlapping promoter. 

Here we analyze how mutual exclusion between tran¬ 
scription factors and RNAP influences the noise level. 
By taking the multi-step transcription initiation explic- 
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itly into account, our study emphasizes that although the 
activity of a promoter may be limited by a single bottle¬ 
neck process, it does matter whether this limiting process 
is early or late in the transcription initiation process. 

II. MODEL 


initiation complex and influence the subsequent RNAP 
activity. This gives the average initiation time interval 
under repressor 1 /r rep ressed = (1 + 1 /K) ■ (1 /r on ) + 1 /rf 
(Fig. [ 2 ]), where the dissociation constant of the repressor 
K = k u /kb quantifies the binding strength of the repres¬ 
sor. The average mRNA number (to) is then given by 

i 


Figure [ 2 ] illustrates the interplay between a simple 
transcriptional repressor, acting solely by promoter oc¬ 
clusion, and the activity of the promoter it regulates. 
The transcription factor binds to the promoter with a 
rate kb when it is free, and unbind with a rate k u . We 
assume the McClure three-step promoter model [Ill, [HI 
for the transcription initiation. The RNAP binds to the 
free promoter with a rate n to form a closed complex. 
Subsequently it can unbind with a rate r_i or form an 
open complex with a rate r 2. The latter step is a non¬ 
equilibrium step, followed by a subsequent elongation ini¬ 
tiation with a rate r 3. 

If there is no transcriptional repression, the total time 
between subsequent promoter initiations can now be ob¬ 
tained by adding together the times for the individual 
steps in the initiation process. This is illustrated in 
Fig -111 where this total time 1 /r is given as the sum of an 
effective on-time 1 /r on = l/ri + (l/ri)-(r_i/r2), and the 
time needed for the subsequent step I/77 = l/r2 + l/r3. 
Noticeably this sum rule incorporates all the three stan¬ 
dard step of the McClure promoter model, with the addi¬ 
tional caveat that the reversible binding step takes some 
additional time because the RNAP may bind and un¬ 
bind several times before the irreversible open complex 
is formed. The total time between two transcription ini¬ 
tiations is accordingly 


1 _ i 1 

r r on 77 


(1) 


where the 0 . 5-1 seconds time interval it take the RNAP 
to move away from the promoter after transcription ini¬ 
tiation for simplicity is included in 1 /r/. Therefore, a 
promoter that is limited by a small elongation initiation 
rate 77 can have an “on-rate” r on which is much hi ghe r 
than its overall initiation rate r = r on ■ r//( r on + 77) | 18 | . 

It should be noted that such multiple sequential steps 
in promoter initiation can reduce the Fano factor below 
one, if each step takes a similar time scale [§). Obviously 
the process will be well-approximated by a single Poisson 
process if one of the steps is the rate limiting, which gives 
the Fano factor one. This is the largest Fano factor that 
the system with the sequential reaction steps can achieve. 

Repressors make the reaction steps to branch out to 
a repressed state, which allows Fano factor to exceed 
one. In Fig 2 , this branching happens around the “/” 
state, with the left branch representing the repressed 
state (“T” state). Noticeably, a repressor that exclusively 
acts through promoter occlusion only interferes with the 
on-rate r on [ 9 j. In other word, when the RNAP is already 
on the promoter, then such a repressor cannot access the 
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where the aspect ratio R = r on /rf characterizes the pro¬ 
moter architecture [ 0 |, and 7 is the mRNA degradation 
rate. 
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FIG. 2. Three-step promoter model of 0 exposed to 
a repressor. The appropriate states is marked T, /,c and 
o, and the figure illustrate how this can be is simplified to 
a process that focuses on the difference between the time 
1/r on of the RNA polymerase association and the time con¬ 
sumed by subsequent steps. In vitro data for LacUV5 is 
ri/r_i = 0.16[RNAP]/nM r 2 = 0.095/sec, r 3 = 2/sec [l3| 
where [RNAP] is free RNA polymerase concentration. 


The described reaction scheme (Fig. [ 2 J provides a 
stochastic mRNA production process. Combined with 
the mRNA degradation at a constant rate 7, the variance 
of mRNA number in the steady state, a 2 = (to 2 ) — (to) 2 , 
can be calculated by using the master equations. We 
performed the calculation for both the full three-step ini¬ 
tiation model and the effective two-step initiation model 
described by the irreversible binding with r on and sub¬ 
sequent elongation with a rate 77 [ 9 |. The detail of the 
derivation is given in appendix [A] We focus on the Fano 
factor v = <t 2 /(to) as a measure of the cell to cell hetero¬ 
geneity, which should be one if the mRNA production is 
a single step Poisson process. 


III. RESULTS 


The Fano factor for the effective two-step initiation 
model with repression, eq. ( 1 A 2 I) . becomes 



{i’onlk b ) - R ■ K • (1 + K) 


[1 + K{ 1 + R )} 2 


( 3 ) 
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FIG. 3. Fano factor as a function of mRNA number (m) with 
r /7 = 15.7, 7 = In(2)/(117sec) ;19] and R = 10, 1, and 0.1. 
Solid lines are for the two-step model (IA2I) , while symbols are 
obtained by three-step model ED with combinations of n, 
r_i, r 2 , and r$ corresponding to the used R and r. Symbols 
with error-bars are the corresponding experimental data from 
Fig. 3 in Jones et al. 0 - Insets show the Fano factor vs. (m) 
where the repressor number fluctuation is taken into account 
by simulating the stochastic production of repressor mRNA 
and protein, which makes kb a stochastic variable. The detail 
of the simulation is given in Appendix [B] (a) Assuming 9.8 
tetramers per cell we set kb = 1/(4.3sec). k u is varied to 
change (m). Inset: Effect of cell-to-cell variation of kb due to 
stochastic repressor production, with 9.8 tetramers per cell 
on average, (b) k u = l/(560sec) from J20j] and kb is varied 
to change (m). Inset: Effect of cell-to-cell variation of kb, 
where the average repressor number is controlled through the 
transcription rate of the repressor mRNA. 


when mRNA degradation rate (typically ~ 1/(3 min) 
[ 2 ll p is much smaller than the repressor binding rate as 
well as the RNAP elongation rate (7 kb and 7 -C 77). 

The importance of the on-rate r on for the cell to cell 
variability becomes evident when we consider the sub¬ 
stantially repressed genes, or genes where the concentra¬ 
tion dependent on-rate of the repressor kb is much higher 
than the off rate k u , i.e., K —> 0. In this case, we have 


[ 9 ]. The difference of eq. © from the simple promoter 
model S0 is that noise can be large for a weak pro¬ 
moter in case its low basal activity is caused by limiting 
later steps in the transcription initiation (e.g. the lac 
promoter). 

Experimentally, the noise is typically measured as a 
function of the average mRNA number (m) i 0, S3- 
The average mRNA number can be controlled by either 
changing the repressor binding strength to the promoter 
(typically by altering the binding site sequence on the 
DNA), or by changing the concentration of the repressor. 
The former corresponds to changing k u at a fixed kb, 
while the latter is the other way around. 

For a fixed repressor concentration (constant kb), from 
© we can express k u as a function of (m). Replacing k u 
in the Fano factor, one get 


v 




(5) 


with m m ax = r/ 7 , where the approximation ignores a 
reduction term in Fano factor, which is small when (to) < 
m m ax, see eq. ©• 

The prefactor is governed by the 7 -corrected associa¬ 
tion rate kl « kb +"f(R+ 1)(1 — Jf n ' > ) ~ kb- This means 
that v decreases monotonically with (m) when the change 
is caused by increased k u (by operator mutations). 

As an illustration we now consider the substantial Fano 
factors that was measured on the Lac system by M- A 
one-step model (R -C 1) would require kb values that 
are much smaller than the overall initiation rate r (~ 
1/(11 sec) for the measured Lac system [Io|,[l9[) to have 
a Fano factor substantially larger than 1, because r on ~ 
r when R -C 1. Indeed :iQ] uses the binding rate for 
one Lac tetramer to be one per 6.3 minutes to fit the 
measured v with a one-step model. However, this rate 
may be too slow given that the association rate of one 
Lac-dimer is estimated to be about 1/3.5 min [2211 . a nd is 
found to bind 5 fold slower than a Lac-tetramer [23l . [ill , 
suggesting an association rate per tetramer of 1/42 sec 
in an E. coli cell. 

The multi-step models can give high Fano factors at 
much higher values of kb . Spassky et al. [Ill measured 
that open complex formation takes l/r 2 ~ 10 sec for 
the lacUV5 promoter in vitro, which combined with r ~ 
llsec suggest that this later step is rate limiting and that 
R 1. Our analysis assuming R = 10 and on-rate of 
a single Lac-tetramer of 1/42 sec gives the Fano factors 
of ~ 4 with ~ 10 tetramers per cell, consistent with the 
experimental data by Jones et al. [Ic| (Fig©i). 

Consider now a given operator (constant k u ) and 
change (to) by regulating the repressor concentration 
(kb). The Fano factor in this case is 


1/ « 



(4) 


V 




( 6 ) 


The increase of v with r on /kb reflects the number of tran- with the 7 -corrected dissociation rate k* ~ + 7 • 

scription initiations between each repressor binding event ( m)/m rnax ~ k u . Eq. © is non-monotonic, with largest 
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v at half maximum expression (m) ~ m max /2 (Fig. ©d). 
The functional dependence of v with (to) in Eq. © does 
not depend on R, but the interpretation of the underlying 
dynamics does. Noticeably, to obtain a given repression 
level rn/rrimax for a promoter with R 1 the repressor 
needs a factor (1 + R) stronger binding than naively ex¬ 
pected. This reflects that the repressor has to act in the 
reduced time where the promoter is not occupied by RNA 
polymerase Q, see Fig. 1. A corollary of this interplay 
is that estimates of repressor binding energies from pro¬ 
moter activities also rely on the non-equilibrium aspects 
of the RNAP-promoter dynamics. 

Finally one may notice that the fit in Fig. 3b underesti¬ 
mates the measured noise level. Part of this is attributed 
to the fact that we use a rate kb that is the same for 
all cells at a given repression level. The Lac-tetramer 
in fact comes at small numbers [25|, even at the highly 
repressed state (to the left of Fig. 3b), and the cell-to- 
cell fluctuations can be substantial. This will add to the 
Fano-factor and in particular increase the variation of 
mRNA for intermediate repression levels, and make the 
maximum Fano-factor to be reached substantially below 
the (to) ~ nimax /2 value predicted by eq. © 

We simulated this effect by considering stochastic pro¬ 
duction and degradation of the repressor (Fig. ©r in¬ 
set). The system is more sensitive to this fluctuation 
for smaller R, because the repressor is more effective 
(see eq. © . In the simulation, we set the parameter so 
that one repressor is produced per one repressor mRNA. 
The effect will be naturally larger if more than four LacI 
monomers are produced per mRNA. Also the effect will 
be stronger when the chromosome copy number is larger 
than one, because fluctuations of repressor act simulta¬ 
neously on each gene copy. 

For Fig. ©i, the fluctuation of repressor number also 
increase the Fano factor (Fig. ©i inset), but a larger R 
shows a larger Fano factor for the same (m). This is 
because to achieve the same (to), k u needs to be smaller 
for larger R , which makes the dynamics more bursty and 
sensitive to the repressor number fluctuation. 


scales of transcriptional repressors in prokaryotes. As 
an added benefit, the formalism propose to use measure¬ 
ments of the Fano factor as a tool to determine the ratio 
of two competing rates (eq. ©). By exposing for exam¬ 
ple a promoter with large r on to different repressors, one 
may compare repressor dynamics. Conversely, by expos¬ 
ing different promoters to the same repressor/operator 
combination, one may quantify their relative on-rates for 
RNAP. To be truly useful, such an experimental design 
should preferentially use a repressor which exhibits min¬ 
imal noise, as one thereby reduce the extrinsic noise. 

Finally, although Fano factors in principle are robust 
to having multiple copies of a given promoter in the E.coli 
cell, then one should be aware that failure in detecting 
all mRNA will make the experimentally measured Fano 
factor systematically smaller than the real one, 

Vmeasured = 1 b P ’ (^ 1) (7) 

where p is the probability for observing a mRNA in the 
cell (detail in appendix O). For instance, the procedures 
based on counting individual spots tend to underestimate 
the number of mRNA molecules 26|; the highest value of 
mRNA per cell reported in ref. [5J, which uses the count¬ 
ing method, is less than 10 , while ref. 6 | that uses the 
total intensity to estimate the mRNA number reports ~ 
50 mRNAs per cell. Thus, if p is say 0.2, then a real burst 
size of v ~ 9 would only be detected as Vmeasured ~ 2.6. 
Therefore a measured Fano factor should be corrected by 
the estimated likelihood for identifying individual mR¬ 
NAs in the cell. 

Using v as a experimental tool to learn about pro¬ 
moter dynamics would further be facilitated by reporter 
mRNAs with relatively large lifetimes (small 7 ). Cen¬ 
tral in such an analysis is to realize that transcriptional 
noise is primarily sensitive to the first steps of the tran¬ 
scription initiation process (Fig. ©, and thereby cell to 
cell variations becomes sensitive to the limiting process 
of individual promoters. 


IV. DISCUSSION 


Appendix A: Derivation of the Fano factor 


The above analyses only apply to repressors that act 
by simple occlusion, and do not affect the post binding 
steps of transcription initiation. In case a transcriptional 
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repressor acts by stalling the isomerization step 
it does not occlude the RNAP binding site and the noise 
should scale with r as suggested by the limit fe). 

In case the transcriptional regulator is an activator, it 
may act through modification of ri, r_i, r 2 or r 3 [lij but 
will not occlude the promoter, and we therefore expect 
the burstiness to be reproduced by considering an overall 
initiation rate modulation as implied in the formalism of 


This short paper aimed to clarify the interplay be¬ 
tween time-scales of transcription initiation, and time- 


We summarize the derivation of the Fano factor for 
the model described in Fig. © In the three-step tran¬ 
scription initiation model, the promoter can be in one 
of the following four states: free (/), RNAP forming a 
closed complex (c), RNAP forming an open complex (o), 
and bound by the transcriptional repressor (T). In this 
model repressor binding does not influence open complex 
formation or the rate of elongation initiation. We denote 
the probability for the promoter to be in the state a and 
having to mRNAs at time t to be P^(t), where a can 
be /, c, o, or T. Assuming that a mRNA is produced 
at the moment the RNAP elongates (this ignores the de¬ 
terministic clearance time), we have the following master 
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equations: 


p m{t) = r 3 P m _i(t) + k u P m {t) 


+7 [(to+ l)P,£ + 1 (f) - 


-1 P m(t) - in + h) P m(t) 
P m(t ) ) 


p m(t) = n p m{t) - (r_i + r 2 )P^(t) 

+7 [(w + l)P^ + 1 (t) - mP^(t)] , 


P° m {t)=r 2 P c m {t)-rzP° m {t) 

+7 [(to + l)P„°+i W - toP£(£)] , 
-Pm(i) = - k u pT{t) 

+7 [(to + 1 )P^+i(t) - mPT(t)\ . 


The probability to have to mRNAs in the system at time 
t irrespective of the promoter/operator state is given by 
p m(t) = Pl(t) + P^{t) + P£(t) + pT{t). The Fano 
factor v = ((to 2 ) — (to) 2 )/(to) was obtained by calcu¬ 
lating (to) = J2m=o mP m and (to 2 ) = I]^ = 0 TO 2 P m in 
the steady state using the generating function method 
The resulting Fano factor for the three-step model 
is given by 


^ 3 —step 


kb 


-K(l + K*) (r 2 a + l2n- + 

V / y P3 rir 2 


7 r 0 


-U 


rlr 2 r 3 / P2P3 


KK* 


K* (R + ^) + ( 1 + (1 + K*) (l + £“)] • [1 + K • (1 + R)] 


(Al) 


with K* = K + 7 /kb and the on-rate r on = ry ■ r 2 /(r_i + 
r 2 ) that is modulated from the rate 77 because the RNAP 
may unbind from the promoter. 

The Fano factor for the effective two-step initiation 
model (RNAP binding and elongation initiation) can also 


be obtained similarly, or by taking r 2 ^ 00 limit of 
eq. ED noting that r on —>• 77 and r/ —> r 3 in this limit. 
The full expression of the Fano factor for the effective 
two-step model is given by 


(r on /k b )-R-K-(l + K*) 

V2 ~ step [1 + K*{1 + R) + ( 7 /r/)( 1 + K*)} ■ [1 + K ■ (1 + R )}' 


(A2) 


Appendix B: Stochastic fluctuation of the repressor 
number 

We assume that the repressor mRNA is transcribed at 
a constant rate a and degraded at a rate T m per mRNA. 
Each mRNA is translated at a rate /3 to produce a re¬ 
pressor and the repressor is degraded at a rate per 
repressor, as parametrized inQ. For simplicity, we as¬ 
sume that a produced repressor corresponds to a LacI 
tetramer. We employ the two-step model in Fig. [2] for 
the promoter dynamics, with making the repressor bind¬ 
ing rate dependent on the number of repressor molecules 
N r as kf, = ko ■ N r , where ko = 1/(42 sec) is the single 
repressor binding rate. For all the simulations, we as¬ 
sumed T m = ln(2)/(120sec), = ln(2)/(40mm), and 


ft = r m , i.e., one repressor tetramer is produced per 
mRNA on average. For Fig. (3K a ) inset, a is set to be 
9.8 • r m ■ T p /P to have 9.8 repressors (tetramers) on av¬ 
erage, while Fig.[3|b) inset, a is changed to control (to). 
The reactions were simulated by Gillespie method |28| . 
and averages were calculated from the data. 


Appendix C: Effect of limited detection on the 
measured Fano factor 

Suppose when we make the observation, each mRNA 
can be observed with a constant probability p. When the 
probability to have to mRNA is P m , then the probability 
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Q(n) to observe n mRNAs is 

Q(n) = V m ~ p n (l—p) m ~ n P m Q(m — n), (Cl) 
^ n\(m — n)\ 

m— 0 ' 

where Q(x) is the Heaviside step function. This gives 

OO 

(n ) = nQ(n) = p(m), 

n —0 

and 


(n 2 ) = p 2 (m(m — 1)) +p(m). 


This results in the measured Fano factor to be 

(n 2 ) - (n) 2 (to 2 ) - (to) 2 

V measured, — / , — P 1 \ T (1 P) 

(n) (to) 

= l+p-(i/-l), (C2) 

where v is the actual value of the Fano factor. 
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